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ABSTRACT 



^ ' We investigate the circular polarization (CP) from Gamma-Ray Burst (GRB) 

CN ; afterglows. We show that a tangled magnetic field cannot generate CP without 

(-\| ' an ordered field because there is always an oppositely directed field, so that no 

■ an ordered field, which carries valuable information on the GRB central engine. 
^ ' By solving the transfer equation of polarized radiation, we find that the CP 
O ■ reaches 1% at radio frequencies and 0.01% at optical for the forward shock, and 

■ 10-1% at radio and 0.1-0.01% at optical for the reverse shock. 

2 ' Subject headings: gamma rays: bursts — gamma rays: theory — polarization 

^ ■ radiation mechanisms: non-thermal — shock waves 

> 

■ 1. INTRODUCTION AND SUMMARY 
, 

Recently a very large linear polarization (LP), ~ 80 ±20%, in the prompt 7-ray emission 
of GRB 021206 was discovered (Coburn & Boggs 2003). The degree of LP was at the the- 
oretical maximum of the synchrotron emission, which implies an uniform ordered magnetic 
field over the visible region (but see also Eichler & Levinson 2003 for the scattering origin of 
LP). Since the causally connected region is smaller than the visible one (Gruzinov & Wax- 
man 1999), an ordered field could be advected from the central engine of the Gamma- Ray 
Burst (GRB) and even drive the GRB explosion (Coburn & Boggs 2003; Lyutikov, Pariev 
& Blandford 2003). 

On the other hand, LP of < 10% (typically a few %) has been detected in the GRB 
afterglows (Covino et al. 2003 and references there in), which is attributed to synchrotron 
emission behind a shock (e.g., Meszaros 2002). In most popular models (Gruzinov 1999; 
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Ghisellini & Lazzati 1999; Sari 1999; Rossi et al. 2002), the magnetic field is generated 
at the shock front and completely tangled (Medvedev & Loeb 1999). LP arises due to the 
geometric asymmetry provided by the afterglow jet observed off-axis if the magnetic fields 
parallel and perpendicular to the jet have different strengths. The 7-ray LP mentioned above 
could be also explained in this model if the jet is very narrow (Waxman 2003; Nakar, Piran 
& Waxman 2003; but see also Granot 2003). 

Thus the present issue is whether an ordered magnetic field exists or not. If an ordered 
field exists in afterglows, its fraction to a tangled field carries valuable information on the 
GRB central engine. In this Letter wc show that observations of the circular polarization 
(CP) could be a useful probe of the ordered field. CP has been detected in AGN jets (Wardle 
et al. 1998; Homan & Wardle 1999; Bower et al. 2002) and X-ray binaries (Fender, et al. 
2000, 2002) in recent years. Theoretically, these observations are explained by a plasma 
effect in synchrotron sources. We apply this theory to the GRBs for the first time. 

There are two main mechanisms to generate CP: intrinsic CP of synchrotron emission 
and Faraday conversion (FC) in sources. FC is a plasma effect which converts LP into CP 
(e.g., Jones & O'Dell 1977a,b). These are treated all together by solving the transfer equation 
of polarized radiation in §2. Then, we show that the tangled field cannot generate CP in §3. 
Next we estimate CP from GRB afterglows in presence of an ordered field together with 
a tangled field in §4. We find that, if the ordered field is comparable to or more than the 
tangled one, the degree of CP is about 1% at radio frequencies and 0.01% at optical for the 
forward shock, and it reaches 10-1% at radio and 0.1-0.01% at optical for the (early) reverse 
shock. The radio CP of the reverse shock remains to be ~ 1% even if the ordered field is 
weak (1% of the tangled one). CP in reverse shock radio emission of GRB 990123 has an 
upper limit of 37% at the 99.9% confidence level (Kulkarni et al. 1999; see also Finkelstein, 
Ipatov & Gnedin 2002). Further observation of CP will be a diagnosis of the ordered field 
and bring clues to the nature of the GRBs. 



2. TRANSFER EQUATION 

We first consider the evolution of the Stokes parameters, /, Q, U, V, (in units of erg 
s~^ cm~^ sr~^ Hz~^) in a homogeneous plasma with a weakly anisotropic dielectric tensor 
(Sazonov & Tsytovich 1968; Sazonov 1969a,b; Jones & O'Dell 1977a; Melrose 1980). It is 
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described by the transfer equation of the polarized radiation, 
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where is the azimuthal projection angle of the magnetic field on the plane perpendicular to 
the line of sight (see Figure 1), K(/,Q,(7,y) are the absorptivity, ^7(/,Q,(7,y) are the emissivity, and 
k'q and Ky are rotativity and convertibility, respectively. These coefficients for a relativistic 
plasma with a power-law energy distribution have been derived by Sazonov (1969a,b) for a 
particular frequency region, Vj^in <^ ^- We have extended the frequency region to <^ ^^mm; 
and both cases are summarized in Appendix A. For the emissivity ?7(/,Q,(7,y), we consider 
only the synchrotron emission. 

There are mainly two ways to generate the circularity V from synchrotron sources. The 
first one is intrinsic emission due to the emissivity r\y (Legg & Westfold 1968). The second 
is FC (e.g., Jones & O'Dell 1977a), which is the conversion of Q and U into V by means 
of Kq in equation (1). If the natural modes of a plasma are nearly circular, the LP vector 
of propagating radiation rotates since the left- and right-circular modes have different phase 
velocities due to birefringence. This effect is well known as Faraday Rotation (FR) (Rybicki 
& Lightman 1979). FC is a similar phenomenon caused by birefringence of a medium. If 
the natural modes are linearly or elliptically polarized, the difference in the phase velocities 
leads to the cyclic conversion between CP and LP. FC becomes effective around the self- 
absorption frequency since \k*qI ki\ ~ 1 for p ~ 2 from equations (Al) and (A3). Note 
that, if the magnetic field is uniform, the synchrotron emission does not generate U in the 
coordinate system with = 0. Therefore FC does not occur without the rotation of Q into 



When we consider the tangled magnetic fields later, the coefficients in equation (1) 
are averaged with respect to the distribution of the magnetic field. This is justified when 
the typical scale over which the Stokes parameters change is much larger than that of the 
field orientation and hence the correlations between the Stokes parameters and the transfer 
coefficients tend to zero (Ruszkowski & Begelman 2002). ^ 



^In this Letter we also neglect the coupling of the natural modes in an inhomogeneous medium (Jones & 
O'Dell 1977b; for its validity see Ruszkowski & Begelman 2002). 
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3. MAGNETIC FIELD CONFIGURATION 

Since wc know very little about the magnetic field configuration in the afterglows of the 
GRB's, we here consider two components, i.e., the ordered field and the axisymmetric tangled 
field. There is a preferred direction in which the fiuid moves. Therefore, the assumption 
of the axisymmetry of the tangled field is quite natural. The ordered field is expected in 
the patchy model (Gruzinov & Waxman 1999), or if the magnetic field advected from the 
central engine prevails or an ordered magnetic field exists in the medium into which the 
shock propagates (Granot & Konigl 2003). The tangled field could be produced by the 
Weibel instability (Medvedev & Loeb 1999) or the turbulence behind the shock. 

We set a coordinate system in the shocked fiuid frame as shown in Figure 1. The z-axis 
is the radial direction, in which the fiuid moves. The ordered field is characterized by the 

strength Bord and the direction {'&ord,y^ord)- The axisymmetric tangled field is described by 
the field strength as a function of i^tan, Bfani'&tan) , and the probability per unit solid angle 
f{^tan)- According to Sari (1999), we adopt Bfani^tan) oc (^^sin^^^to,, + cos^ t^to,)"^/^ ^nd 
fi'&tan) oc B^^^i^tan)- If ^ 1) {Bf) ^ {B\) , and vice versa, where B\\ and -B_l are the 
tangled magnetic field components parallel and perpendicular to the z-axis, respectively. We 
parametrize the ratio of the ordered field to the tangled one by C = B'^^/{{Bp + (-B^)). The 
total strength B"^ — B"^^ + (Sy) + {BD is determined by the afterglow model in the next 
section. 

Given a magnetic field geometry and the plasma parameters from the afterglow model, 
we can take an angular average of the coefficients and solve the transfer equation (1) to 
obtain the Stokes parameters. Although the shock structure may be important (e.g., loka 
2003), we neglect it as a first step. We define the observed frequency as the frequency in the 
fiuid frame multiplied by the Lorentz factor of the fiuid in the lab frame 7, whereas ratios 
of the Stokes parameters are Lorentz invariant. 

If C = (no ordered field), we can immediately find that (77Qsin20), (rjv), (KQsin20), 
{kq sin 20) , (kv) and (Ky) vanish as a result of the axisymmetry and refiection symmetry 
about the xy plane. Therefore the tangled field alone cannot generate CP. Intuitively this is 
because there is always an oppositely directed pair of magnetic fields, so that no handedness 
exists. 



4. APPLICATIONS TO GRB AFTERGLOWS 

First we consider the forward external shock with energy E propagating into a constant 
surrounding density n. According to the standard afterglow model (Sari, Piran & Narayan 
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1998), the Lorentz factor of the shocked fluid and the radius of the shock evolve as 7 = 

respectively, where t is the observer 
time. The shell thickness in the shocked fluid frame can be estimated by i?/7, which we use 
as the path length of the transfer equation (1). We assume that electrons are accelerated in 
the shock to a power law distribution of Lorentz factor 7e, iV(7e)(i7e oc %^dje for 7e > 7mm, 
where = , N{'fe)d'ye = 47n is the electron number density in the shocked fluid frame 
and p > 2. The minimum electron energy and the magnetic field strength in the shocked 
fluid frame are given by 7^j„ = eg [(p — 2)/(p — 1)] {mp/me)^ and B — {32Tr'mpeBny/^jc, 
respectively, where the parameters ee and eg are fractions of shock energy that go into the 
electrons and the magnetic energy, respectively. Thus we can calculate 7, 7mm, B, R/^ and 
A'^e as a function of t for given E, n, p, and e^. Wc adopt E = 10^^ erg, n = 1 proton cm~^, 
p = 2.2, Ce = 0.1 and — 0.01 (Panaitescu & Kumar 2002). For simplicity, we temporarily 
neglect the jet effects. 

The typical synchrotron frequency is i^m = ei?7^j„7/27rmeC = 145eg _;^e^^^^2-^52^^iday 
GHz where the convention Q = lO^Qx is used except for tiday — t/1 day. For v > and 
p < Urn, we use equations (Al) and (A3), respectively. By using equation (A3), we can 
estimate the self-absorption frequency from KjR/'y = 1 as = 41e7i^e]/^2-^52^^^''^ GHz, 
where we put sin^^ = 1. The flux is given by F,y oc t^^'^u^ with (3 = 2 and 1/3 for i/ < 
and Va < V < Vm-i respectively (Sari, Piran & Narayan 1998). We can neglect the electron 
cooling for our interests. 

Figure 2 illustrates CP from a forward shock. The degree of CP is about 1% at the 
self-absorption frequency Va and 0.01% at optical for C ^ 1- We have found that the de- 
pendence of CP on t and ^ is weak, and the dependence on C, i^ord, ^ord and % is roughly 
oc •\/C/(l -I- C) cos 9ord (see Figure 1 for the relation between 9ord, '&ordi Vord and %). The CP 
is mainly due to intrinsic emission. 

Next we consider the reverse shock propagating the ejected shell itself (Meszaros & Rees 
1997; Sari & Piran 1999a). Let 70 and T be an initial Lorentz factor of the shell and the 
burst duration, respectively. Then the Lorentz factor at the shock crossing time is given by 
7x = min(7o,7c), where 7c = (3£'/327rnmpC^T'^)^/^ is a critical Lorentz factor (Kobayashi 
2000; Kobayashi & Zhang 2003; Zhang, Kobayashi & Meszaros 2003; Sari & Piran 1999a,b). 
The shock crossing time is given by tx = {lo/lx-Y^^T . At t = tx the minimum electron 
energy, the magnetic field strength and the electron number density in the shocked fluid 
frame are given by j„iin = [{p - 2)/ {p - 1)] (mp/me)7o/7x, B = (327rmpeBn)^/^7xC and 
Ne = 47xn/7o, respectively. After the shock crossing, these quantities approximately evolve 
as 7 oc r^/^^ 7mm OC r^^/^^ B oc r^3/2^ and Ne oc r^^^^^, respectively (Kobayashi 2000; 
Kobayashi & Sari 2000). The shell thickness in the shocked fluid frame is R/^ where R = 
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{nEt/inmpncy/'^. We adopt 70 = 100 and T = 100 s. 

The typical synchrotron frequency is = 0.035 el_i £5^2^^'^^ hJay^^ 2^73/48 ^2^^ 
max[l, (7o/7c)^^^^^] GHz. Using equation (Al) we can estimate the self-absorption fre- 
quency from kjR/^ = 1 as = 20e^;!7^)/^+^) eg+f'^^"^ n(^+^y'(^+^) hZ?^'''^^''^-''^ 
y^(73p+i46)/48(p+4) £;V2(f+4) max[l, (70/7^) GHz, where we put sin^ = 1. 
The flux peaks at ~ We can neglect the electron cooling for our interests. 

Figure 3 shows CP from a reverse shock. For C ^ 1) the degree of CP reaches 10-1% 
at radio frequencies, and 0.1-0.01% at optical. Even if C ~ 10~^, CP remains to be ~ 1% 
at the self-absorption frequency Ua- This is because of the FC. The intrinsic CP decreases 
as C decreases, but V is generated from Q and U due to the FC. The dependence of CP 
on ^ is weak. If FC is not effective, the dependence on (, i^ord, Vord ^-nd x is roughly 
-\/C/(l + C) cos Oord, while if FC is effective, it is not so simple. 

So far we have calculated CP at one point on the afterglow image on the sky. The 

integration over the entire emitting region may suppress the observed CP as in the case 
of LP (Sari 1999; Ghisellini & Lazzati 1999). Let us estimate the suppression factor when 
we observe a jet with an opening angle ©q from a viewing angle G„. We consider the 
ordered magnetic field defined by ("i^ordj Vord) in the coordinate with x-axis being the direc- 
tion from the jet center to the line of sight on the sky and with 2;-axis being the direction 
in which the ejecta moves (see Figure 1). We assume that the viewable region is a uni- 
form disk of angular extent I/7 centered around the line of sight to the observer, since 
the afterglow image is rather homogeneous for v < v^^ (Granot, Piran & Sari 1999a,b) 
or after the jet break (loka &: Nakamura 2001). Then, if we assume V/I ex cos 9ord and 
take the coordinate origin at the line of sight on the sky, the suppression factor is given by 
J dQd^ [cos i9ord cos xe + sin xe sin 19 ordCos(ipord — ^)] / J dQd^, where sinxe = 276/(1 + 
7^0^) and 0^ + 20^,0 cos $ + 0^ < 0g. Interestingly, no cancellation takes place when 
''^ord = or TT. Even when the cancellation occurs, e.g, in the case of {}ord = 7'"/2, some 
amount of CP remains if the visible region has an asymmetry due to the jet geometry as in 
the case of LP. 

We are grateful F. Takahara for useful comments. This work was supported in part by 
Grant-in- Aid for Scientific Research of the Japanese Ministry of Education, Culture, Sports, 
Science and Technology, N0.OO66O (KI). 
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A. Transfer coefficients 

We summarize the transfer coefficients in equation (1) at an angle 9 to the magnetic field 
B (Sazonov 1969a,b; Sazonov & Tsytovich 1968; Melrose 1980). In this section we measure 
all quantities in the shocked fiuid frame. We assume that the electron number density in the 
interval of the Lorentz factor (i7e is power-law dN^ — N{'^e)dle — -^eT^^c^Te for Imin < 7e 
with p > 2. Then, at frequencies Vmin = 7mm ^-L ^ ^ coefficients are given by 

= /^,K^(z/K)-°-5/^ = «:y = «:^«:±(i/K)-"-'cot^, 

K*y = K*jK^{u/u^)-\ln^rmnh;;^it^'^ cote, (Al) 

where i'± = |e|Ssin^/27rmeC, 'r]± = {e^ /c)Nei'±, k±_ = {e^/mec)Ne/i'±, a = {p — l)/2, and 
3«+'/^ ^/a n\^fa 1\ Q 3«+"' ^fa fa 1 
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The above notations are the same as Jones & O'Dell (1977a; But k*^ is different. The 
appropriate integration of equation (9) in Sazonov (1969a) gives the above value). For 
the frequency region, v ^ Vmin-, the analogous representations have not been derived yet. 
Such a frequency region becomes important in the application to the forward shock of GRB 
afterglow. In this case, we obtain the following expressions: 



I I \l/3 -2a-2/3 1 V -2a-l j. n 



/ I \-5/3 -2Q-5/3 1 y / / \-2 -2(a+l) , 

Ki = K^t^l.[ylVl) ' Irrnn . I^Q = l^t^I, /^V = «± (^^/^^x) 7™„ COt I 



4 = <^'^±(^/^±)-'/'7;.r'^', /.^ = <^/.±(^/^±)-^(ln7m.n)7;.^^'^cot^, (A3) 
where 

3'/" (2\ 3V=« + 3/2 m 

2(a+l/3) Vsy' ° 3(a + l/2)' ° 2a + hl6 
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Note that the intrinsic LP {rjQ/rji) in this frequency region is not (p + l)/(p + 7/3) 70% 
but 50%. In GRB 021206, considerable photons are below the break energy, i.e., v <^ v^nin 
(Coburn & Boggs 2003), so that even an ordered field might not explain the observed LP 
fa 80 ± 20%. 
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Fig. 1. — Magnetic field orientation in the shocked fluid frame, {left panel) The fluid flow 
is along the z axis. We introduce the ordered and tangled magnetic field. The former is 
defined by (J) ord, V ord) , while the latter is stochastic and distributes as a function of {'dtam 
^tan)- The observer is in the direction of the axis z\ which is specified by the angle % on the 
xz plane, {right panel) It is convienient to take a new coordinate in order to deal with the 
radiation transfer. In this coordinate, the magnetic field direction is specified by {9, (p). 
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Fig. 2. — Circular polarization (= from a forward shock is shown as a function of 

the observed frequency u ior t — 1 and 0.1 day. The self-absorption frequencies Ua are 
marked by black points. The ratio of the ordered field to the tangled one is defined by 
C — ^ord/i{^\\) + (-^i))- The ordered field is directed to 'i9ord = 0. The perpendicular 
tangled component dominates the parallel one, i.e., ^ = 0.1. The observer is in the direction 
X = 7r/4 (see Figure 1). For this configuration, V < 0. We adopt E = 10^^ erg, n = 1 proton 
cm-3, p = 2.2, ee = 0.1 and sb = 0.01. 
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Fig. 3. — Circular polarization (= |V^|//) from a reverse shock is shown as a function of the 
observed frequency z/ for t = 1, 0.1 and 0.01 day. For V < {V > 0) we use solid (dashed) 
lines. The self- absorption frequencies Ua are marked by black points. The dotted lines are 
calculated by putting Kq = Q (no conversion) . The ratio of the ordered field to the tangled 
one is defined by C = -^ord/ + i^l))- We adopt E = 10^^ erg, n = 1 proton cm~^, 
p = 2.2, Se = 0.1, SB = 0.01, 7o = 100, T = 100 s, ^ord = 0, ^ = 0.1 and x = 7r/4- 



